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We propose the integrable A-dimensional Calogero-Coulomb-Stark and two-center Calogero- 
Coulomb systems and construct their constants of motion via the Dunkl operators. Their 
Schrodinger equations decouple in parabolic and elliptic coordinates into the set of three differ¬ 
ential equations like for the Coulomb-Stark and two-center Coulomb problems. The Calogero term 
preserves the energy levels, but changes their degrees of degeneracy. 


I. INTRODUCTION 

One of the most important features of the Coulomb problem is its maximal superintegrability caused by the 
conservation of the Runge-Lenz vector. As a consequence, its Hamiltonian admits separation of variables in several 
coordinate systems, and hence, a few types of integrable perturbations with separation of variables. The textbook 
examples of such perturbations are (see, for example, Ref. 1): 

1. The Coulomb problem in constant electric field (the Coulomb-Stark problem) with the Stark potential 5V{r) = 
F ■ r, which admits a separation of variables in parabolic coordinates. 

2. The two-center Coulomb problem, which admits a separation of variables in elliptic coordinates. 

These systems remain integrable despite the partially broken rotational symmetry along the highlighted direction. 
The latter coincides with the direction of the electric field F for the Coulomb-Stark problem, and with the vector a 
connecting two Coulomb charges for the two-center Coulomb problems. As a result, the orthogonal angular momentum 
components and the modified longitudinal component of the Runge-Lenz vector are preserved. 

However, the aforementioned systems are not solvable exactly. In the Coulomb-Stark problem one can get an¬ 
alytically only the perturbative spectrum, while in the two-center Coulomb system, the energy spectrum can be 
constructed only numerically, except for some special cases [2]. Nevertheless, the separation of variables is crucial in 
the study of these systems. 

Recently, together with O. Lechtenfeld, we have observed that the TV-dimensional Coulomb problem deformed by 
the rational Calogero potential [3] (and by its generalization associated with arbitrary root system [4]) remains super- 
integrable and has the same energy spectrum as the original Coulomb model [5]. The Calogero-Coulomb Hamiltonian 
is given by the expression 
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In Ref. [7] we have revealed an explicit expression for the analog of the Runge-Lenz vector in this system. We also 
observed that being formulated in terms of the Dunkl operators [8], the conserved quantities and their algebra look 
pretty similar to those in the conventional Coulomb problem. This led us to the claim that most of the properties of 
the conventional Coulomb problem and its integrable perturbations have their Calogero-Coulomb counterparts. 

Thus, one may ask the following: 

Are the high-dimensional Coulomb-Stark and two-center Coulomb problems with the Calogero potential (we will 
refer them as Calogero-Coulomb-Stark and two-eenter Calogero-Coulomb problems) integrable systems? If they are, 
do they admit separation of variables, at least partial? 

In this paper we give positive answer to both questions: 
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i) The proper choice of the direction of the vectors F and a makes the related systems integrable. The corresponding 
Hamiltonians are given by the expressions 
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The constants of motion are the Dunkl-operator deformations of the corresponding integrals in the underlying Coulomb 
systems with the subsequent symmetrization over the coordinates. 

Note that the integrability of the Coulomb-Stark and two-center Coulomb systems (even without the Calogero 
interaction, even at the classical level) does not follow from the existence of the Runge-Lenz vector or the properties 
of the symmetry algebra of the usual Coulomb problem. The preservation of the constructed constants of motion in 
the systems (1.2) and (1.3) will be verified independently, based on the properties of the Dunkl operators. 

ii) After the transition to the Jacobi (center-of-mass) coordinates the above listed systems admit a complete sep¬ 
aration of variables in parabolic/elliptic coordinates for N = 2,3 and a partial separation for iV > 3. In fact, the 
Schrodinger equation decouples into three parts, only one of which depends on the Calogero inverse-square term. The 
latter can be treated as a deformation of the Schrodinger equation for the SO{N — 1) angular momentum, usually 
referred to as an angular Calogero Hamiltonian [9-12]. 

The paper is organized as follows: 

In Section 2 we describe the Calogero-Coulomb problem and its symmetry algebra, which contains the deformations 
of the angular momentum and Runge-Lenz vector by means of the Dunkl operator [5, 7]. Then we perform the 
orthogonal transformation to the Jacobi coordinates and show that the transition to the spherical coordinates factorizes 
original Schrodinger equation into three parts. The first two are just the radial and angular Schrodinger equations 
of the standard Coulomb problem. The third part produces the Schrodinger equation of the angular Calogero model 
with the excluded center of mass, which has been solved in Ref. 10. 

In Section 3 we introduce the Calogero-Coulomb-Stark problem and present explicit expressions for its constants 
of motion. They include the components of the Dunkl angular momentum, which are orthogonal to the direction of 
the electric field F and the modified component of the Runge-Lenz vector along the field direction. We go on to show 
that after the transition to the Jacobi coordinates, the system possesses a partial separation of variables in parabolic 
coordinates as in spherical coordinates. Namely, the Schrodinger equation decouples into three parts, where the first 
two coincide with the respective Schrodinger equations of the standard Coulomb problem in parabolic coordinates, 
while the third one produces the Schrodinger equation of the angular Calogero model with the excluded center of 
mass [10]. 

In Section 4 we introduce the two-center Calogero-Coulomb problem and perform the similar study for it. Again, 
the components of the Dunkl angular momentum, orthogonal to the symmetry axis of the system given now by the 
vector a, are preserved. The last integral is a deformation of similar construction applied for the underlying two- 
center Coulomb system [13]. The latter is based on the peculiar mixture of the longitudinal components of the two 
Runge-Lenz vectors, corresponding to each charge. After transition to the Jacobi coordinates, the system acquires a 
partial separation of variables in elliptic coordinates. 


II. CALOGERO-COULOMB PROBLEM 

A. Coulomb problem 

The A-dimensional Coulomb problem is defined by the Hamiltonian 
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It possesses the maximum number (2N — 1) of functionally independent constants of motion. First, due to the 
rotational invariance, the system possesses conserving angular momentum tensor 

L^j — XiPj XjPi- ( 2 . 2 ) 
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Its components satisfy the standard commutation rules for the SO{N) generators: 

\^ij: (2.3) 

The angular momentum provides the system with iV — 1 commuting quadratic (in momentum) constants of motion 

Ll= Llj with 2<k<N, L% = L'^. (2.4) 

They correspond to the quadratic Casimir elements, constricted from the generators of the sequence of the natural 
embeddings 


SO{2) c • • • C SO{k) c • • • C SO{N). 

Together with the Hamiltonian (2.1), these constants of motion form the complete set of the Liouville integrals. 

Apart from the angular momentum, the A^-dimensional Coulomb motion preserves the Runge-Lenz vector defined 
by expression 


= (2.5) 
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where the curly braces mean the anticommutator: {a, b} = ab + ba. 

The angular momentum tensor and the Runge-Lenz vector define the complete set of constants of motion for the 
Coulomb system. 


B. Calogero-Coulomb Hamiltonian 

The Calogero-Coulomb problem introduced in Ref. 6 contains the additional inverse-square interaction term a la 
Calogero, as already mentioned in the Introduction. It is a mixture of the A^-particle rational Calogero model [see 
Refs. 4 and 14 for review] and of the TV-dimensional Coulomb system. In vector notations the Hamiltonian (1.1) can 
be rewritten as follows 
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It inherits most of the properties of the original Coulomb system (2.1) and possesses hidden symmetries given by 
an analog of Runge-Lenz vector [7]. It is convenient to describe this system by means of the Dnnkl operators which 
makes transparent its analogy with the Coulomb problem. Namely, instead of the 'H-y, let us consider the more general 
Hamiltonian 
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where the momentum modified by the Dunkl exchange operators [8] is used, 

TT = -iV, y,=di-Y — z — 


(2.7) 


(2.8a) 


The operator Sij in Eqs. (2.7) and (2.8a) interchanges the T-th coordinate with the j-th coordinate, 

SijN'(. ..,Xi,...,Xj,...) = ^'(. ..,Xj,...,Xi,.. .). (2.8b) 

On the symmetric wavefunctions the generalized Hamiltonian "H®®” reduces to the Calogero-Coulomb Hamiltonian 

( 2 . 6 ). 

The Dunkl operators commute mutually like ordinary partial derivatives. However, their commutations with coor¬ 
dinates are nontrivial deformations of the Heisenberg algebra relations [8] , 

[xi,TTj] = 0, [xi.Xj] = -iSij. (2.9a) 
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Here the operators Sij are related with the permutations in the following way: 


-gSij for i ^ j, 

^ + 9Y.k^^s^k fori=j. 


(2.9b) 


Using the above relations, it can be verified, that the operator-valued tensors [Sij,Xk\ and [Sij,TTk\ are symmetric in 
all three indexes [9], 


It must also be mentioned that the symmetries of the two-dimensional Calogero-Coulomb problem associated with 
the dihedral group D 2 have been recently studied in Ref. [15] in terms of the Dunkl operators. 


C. Constants of motion for 7^®®" 

Let us define the deformed angular momentum operator via the Dunkl momentum [9, 16, 17]: 

Lij = XiTTj — XjTTi. (2-11) 

Its components satisfy the following commutation relations [9], 


“b iLjlSki '^LilSkj '^LjkSli- (2-12) 

The Dunkl angular momentum inherits many properties of the ordinary angular momentum (2.2). In particular, the 
equations 


= 27ri, '^{Sik.Xk} = 2xi, (2.13) 
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which follow from the relations (2.9a) and (2.9b), imply that the operators Lij commute with the Dunkl momentum 
square and with the radius vector length [9]: 


[7r2,Ly] = [r,Ly] =0. (2.14) 

Hence, the Dunkl angular momentum is a constant of motion of the generalized Calogero-Coulomb Hamiltonian, [5]: 

[•Hf",Ly]=0. (2.15) 

The deformation of the Runge-Lenz vector (2.5) is given by the expression [7] 
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It contains the permutation-group invariant element which vanishes in the absence of the Calogero term [g = 0) 

s = J2s.r- [S,S.j] = 0. (2.17) 

i<j 


The deformed Runge-Lenz vector is also a constant of motion: 


(2.18) 


The proof is given in Ref. 7, and we repeat it here for the completeness. 

According to (2.14), the commutator of the Hamiltonian with the first term from the rhs of Eq. (2.16) is proportional 
to 


[Efi*i.v)U] = = {7.-} - . 


(2.19) 
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and its commutator with the last term is proportional to 
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Combining the relations (2.19) and (2.20) and using the definition (2.17), we obtain: 
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In the last equation we used the identity 


'y — [S', TTi]. 
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Finally, Eq. (2.21) completes the proof of the conservation condition (2.18). 


( 2 . 20 ) 
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The deformed Runge-Lenz vector can be written directly through the coordinate and Dunkl momentum by avoiding 
the invariant S: 

A = X ~ ~ - 2 —^ (2.22) 

Here Pr is the ordinary radial momentum canonically conjugate to the radial coordinate: [pr,?*] = —i. It turns into 
its 3 = 0 analog under the replacement tt —>■ p. 

The relation (2.22) can be checked by inserting Eq. (2.11) into Eq. (2.16) and applying the identity 

TT ■ X + X ■ TT = 2{x ■ p) — iN = 2rpr — iN. (2.23) 

This identity is a direct consequence of the deformed canonical commutation relations (2.9a). 

The commutation relations among the Dunkl angular momentum components (2.12) are supplemented by those 
which incorporate the deformed Runge-Lenz momentum and Hamiltonian [7]: 

[H„ Lu] = lAuSu - lAiSu. [A*, A,] = -2inf’^L,,. (2.24) 

In the absence of the inverse-square interaction, the exchange operator is reduced to the Kronecker delta, and we 
get the familiar algebra, formed by the symmetry generators of the Coulomb problem. On a negative constant-energy 
surface = if < 0 it reduces to the orthogonal algebra so{N + 1), and on the positive one, = if > 0, to the 
so{N, 1). 


D. Constants of motion for 

Let us remember that the Calogeo-Colomb problem can be obtained by the restriction of the generalized Hamilto¬ 
nian (2.7) to symmetric wavefunctions (2.8b). Therefore, its constants of motion can be constructed by taking the 
symmetric polynomials depending on the components of the Dunkl angular momentum and Runge-Lenz vector [5, 7]: 

L2k = Y.Ljl (2.25) 

i<j 

Au = Y.AI (2.26) 

i 

The expressions above demonstrate that the Calogero-Coulomb problem is a superintegrable system, like the pure 
Calogero [18] and Coulomb models. Among these integrals only the two are quadratic on momenta. The first one is 
the square of the Dunkl angular momentum. 


£2 = 21 +S{S-N+ 2 ), 


where I is the angular part of the Calogero (Calogero-Coulomb) Hamiltonian [9]: 
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It coincides with the Casimir element of the Dunkl angular momentum algebra and is also a conserved quantity: 




(2.29) 


So it can be used instead of £2 as a constant of motion. 

The second quadratic integral is given by the sum of the components of the deformed Runge-Lenz vector defined 
in Ref. [5]. For convenience, we take a rescaled expression, 

= V) - ™ = »^o(2-Hr" + 2 ) - (rp. + (2.30) 

It depends on the normalized center-of-mass coordinate and momentum 

1 iV ^ N 

= Po = ^Y.P^- (2.31) 

The existence of two second-order integrals suggests that the Calogero-Coulomb problem admits a separation of 
variables for dimension N < S. 


E. Jacobi coordinates 


The Jacobi coordinates [11, 19, 20] guarantee the separations of the center-of-mass (2.31) from the relative motion, 
keeping the kinetic term’s shape. They are defined in the following way: 


yo — xq — -^={xi -\ --I- xn), 
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(xi -I- \-Xk- kxk+i), I <k < N - 1. 


(2.32a) 


The first coordinate describes the center of mass, while the other coordinates characterize the relative motion. It can 
be verified that the above transformation is orthogonal: 


N N-1 

yk = '^ Okixu = X! ^ikyi- 

1=1 1=0 


(2.32b) 


In order to express the Hamiltonian and Dunkl momentum in terms of the Jacobi coordinates, we note that the 
permutation operators are just orthogonal reflections across the hyperplanes o: • a: = 0 in the TV-dimensional coordinate 
space: 


2 ( 0 . ■ x) 

SijX = SaX := X - a with a = otij = — e^, 


(2.33) 


where form the standard basis in : x = The operators Sij are expressed in a similar manner. Note 

that this representation admits the generalization of the exchange operators (2.8a) to an arbitrary finite reflection 
group known also as a Coxeter group [8] . The vector a is called a root of the Coxeter system. A particular choice of 
roots in Eq. (2.33) defines the Ajv-i Coxeter system. They separate the center-of mass coordinate from the relative 
motion. 

We rewrite the reflections (2.33) in the Jacobi coordinates using the vector notations: 


SijV = Spy = y- 


2(/3 ■ y) 

/3-/3 


(2.34a) 


Here the vectors /3 are images of the corresponding root vectors a under the Jacobi transformation (2.32b): 


Pk = ^ Okk'Oik', 

k' 


Pij,k — Oki Okj • 


(2.34b) 


In fact, the reflection (2.34a) is restricted to the relative 


Their explicit form (2.32a) ensures the relation 0 = 0. 
coordinates (see Eq. (2.40) below). 
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In Jacobi coordinates the inverse-square Calogero potential in the Hamiltonian (2.6) reads 






2 ■ 


(2.35) 


Evidently, the Dunkl momenta along the new coordinate directions are related to the usual Dunkl momenta TTa,. = tt^ 
by the vector transformation 
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— '^Vi — 


(2.36) 


i' = l 


Its zeroth component describes the center-of-mass momentum (2.31): = po = Po- In order to ensure the invariance 

of the defining relations (2.9a), the operators Sij must behave as a second-order tensor in an orthogonal transformation: 


N 

— 0; — 'f'Sij — Z ^ ^ Oii'Ojjf Sif j/. (2.37) 

In Jacobi coordinates, the Cherednik algebra decouples into two commuting parts describing the center-of-mass 
(?/o = xq , Po) Slid relative (y^, Tfj with 1 < i — 1) motions, respectively. Indeed, exploiting the defining relations 
(2.9b) and the orthogonality of the Jacobi transformation (2.32b), we obtain 

AT N N 

Sio = -y= ^ Si'ji ^ Oiii = ^ Oii'Ooii = SiQ. (2.38) 

j'=i i'=i i'=i 


It implies 


[^i,yo] = [^i,a;o] = [po,yi] = -iSio- (2.39) 

Note that the constraint (2.38) is equivalent in old coordinates to the condition Sij = 1, which follows from the 
defining relation (2.9b). 

Below we use the following vector notations for the description of the relative coordinates and momenta: 

y={yi,---,yN-i), fr = (Ifi,.. .ttat.i). (2.40) 

Note that the operators do not permute the coordinates any more. However, the tensors [Sij,yk] and [5^,71^] still 
remain symmetric so that the relations (2.10) are preserved: 

[Sik,yj] = [Sjk,yi], [Sik,^j] = [Sjk,T^i]- (2.41) 

This property persists for any Coxeter system [9]. 

Of course, similarly to the operators Sij , the Dunkl angular momentum transforms as a second-order tensor under 
an orthogonal transformation: 

N 

Lij = Ly-y^ = 'y ( OurOjj'Liijt. (2.42a) 


The rotated components obey the commutation relations (2.12) with the permutation generators Sij replaced by their 
shifted counterparts Sij. They acquire the standard form (2.11) in the Jacobi coordinates and momenta, 


Lij = PiTfj — yj^j with 0 < i, j < iV — 1, (2.42b) 

Loi = Xqtt^ - pjpo. (2.42c) 

Using (2.38), we can rewrite them in the form which splits up the center-of-mass and relative degrees of freedom: 

[Lij ^ Lki\ — '^LikSij -j- iLijiSki ‘^LiiSkj zLjkSii^ (2.43a) 

\_Lij^ Tq/] — zI/QjSii tljQiSij^ (2.43b) 

[Loj, Lqi] = iLji with ^ < i, j,k,l < N — 1. (2.43c) 


So the generators (2.42b) form the deformed analog of the so{N — 1) angular momentum algebra consisting of the 
relative degrees of freedom. 

Usually, in the Coulomb problem, the separation of variables is performed in spherical coordinates. However, for 
the (partial) separation of variables in spherical coordinates, we should first perform the orthogonal transformation 
to the Jacobi coordinates in order to extract the center of mass [11]. 


7 




F. Separation of variables in spherical coordinates 


Let us extract now the coordinates of the center of mass from the Calogero-Coulomb Hamiltonian (2.6), 


^7-y- 




no 


with 


= V^- 


(2.44) 


The last term is just the Calogero Hamiltonian with the excluded center of mass, which depends on the iV — 1 relative 
degrees of freedom. In the Jacobi coordinates, it reads: 


Ho 


y- g(g-1) 
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(2.45) 


For the simplest N = 2 case, the relative Calogero model is reduced to a one-dimensional conformal mechanical 
system with the Hamiltonian 


Ho = 


E. 
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g(g-1) 

22/2 


In the TV = 3 case, it takes more complicated form: 

^ Pi +pI g(g- 1) 2g(g- 1) 

2 + 22/2 ^ (732/2+ gi)2 


2 g(g- 1) 
(73^2 - yi)^ 


(2.46a) 


(2.46b) 


In the general case, due to the conformal symmetry, the relative Hamiltonian in (2.44) can be spit further into its 
radial and angular parts. 


2 y2 


(2.47) 


where (p,, are the angular coordinates, parametrizing the (TV — 2)-dimensional sphere, and I is the Hamiltonian of 
the angular Calogero model with the excluded center of mass (relative angular Calogero model). The latter system 
has been investigated in detail in the series of papers [10-12]. For the TV = 3 case it is reduced to a one-dimensional 
exactly solvable system on circle, introduced by Jacobi in the middle of the 19th century [21]: 


9g(g-1) 

2 2 cos2 


(2.48) 


For higher dimensions, the relative angular Calogero Hamiltonian has more complicate form (see Refs. 11 and 22). 

Now we are ready to transition to spherical coordinates. The azimuthal angle 9 is chosen to measure the projection 
of the position vector to the relative coordinate’s hyperplane: 


Xo = rcos9, y = rsinO h{ipP). (2.49) 

Here the (TV — l)-dimensional unit vector n is parametrized by the remaining angles ip,, with 1<6<TV — 2. In these 
coordinates the Hamiltonian (2.28) looks as follows: 


n^ = 


1 9 ^N-i d 

2j,Af-l Qj. Qj. 


I{9,ip,,de,d^J 
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The angular Calogero Hamiltonian I is related with the relative angular Hamiltonian I by the expression 


I{e,(p„de,d^J = 


n{<PL,d^J 
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(2.51) 


Both quantities I and I are constants of motion. They separate the radial r and azimuthal 6 coordinates, each 
being symmetric in the primary coordinates Xi, from the relative angular coordinates (/7t. Let us choose the factorized 
wavefunction: 


'i!{r,e,^p,) = i?(r)4>(6»)i/>(v?J. 


(2.52) 
















It decouples the Schrddinger equation = E'^ into the following system of differential equations: 
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g(g + JV_3) 2 )^ ^ 0 ^ 

Sin 0 J 

(2.53b) 

{2i{+,,d^^) - q{q + N -3)^ ’ipqi+i) = 0. 

(2.53c) 


The last equation describes the spectrum and eigenstates of I, which were recently investigated in detail in Ref. 10. 
In particular, the spectrum is determined by the numbers 

~ ^ + Sh + ... + NIn with Z, =0,1,2 ,.... (2.54) 

For integer values of the coupling g, the angular energy spectrum is that of a free particle with angular momentum q 
on the (N — 2)-dimensional sphere, but it has a signihcantly lower degeneracy due to the restriction to the symmetric 
wavefunctions [10, 22]. 

The first two equations above are similar to those for the A^-dimensional Coulomb problem. The solutions of 
both equations can be found, for instance, in Ref. 23. The second equation (2.53b) determines the spectrum and 
eigenfunction of the angular Hamiltonian I. For integer coupling g, its spectrum coincides with that of a free particle 
with angular momentum I on the {N — l)-dimensional sphere [10]. The first equation (2.53a) is related to the radial 
Schrddinger equation for the A'^-dimensional Coulomb problem. 

For the positive integer values of the coupling constant, g = 1,2,..., the solutions of Eq. (2.53b) exist, provided 
the quantum number I takes the values 


I = q + li with /i = 0,1,2,.... (2.55) 

so that I > ^N{N — l)g. 

Similarly, the well-defined solution of the radial equation (2.53a) exists for the energy spectrum depending on the 
principal quantum number n: 

En = - 7T, n = Ur + I + I, nr = 0,1,2, - (2.56) 

2{n+^)^ 

Hence, the Calogero-Coulomb problem remains degenerate with respect to q. The energy spectrum of the system 
is the same as in the A'^-dimensional Coulomb problem, but with a lower degeneracy due to the restriction to the 
symmetric wavefunctions on the initial coordinates. 


III. COULOMB-CALOGERO-STARK PROBLEM 
A. Coulomb-Stark problem 

It is well known that the TV-dimensional Coulomb problem in constant uniform electric held F remains integrable 
and admits a separation of variables in parabolic coordinates. Its Hamiltonian is dehned by the expression 

Kf = ^--+F-^- ( 3 - 1 ) 

' 2 r 

The perturbative spectrum of this system can be found in Refs. [1, 23]. 

The Stark interaction breaks the SO{N) rotation symmetry of the Coulomb problem down to the SO{N — 1) 
symmetry. The corresponding constants of motion are given by the components of the angular momentum (2.2) 
which are orthogonal to the direction of the electric held, 

L^j = Lij -f TijTikEif^., (3.2) 

Here n is the unit vector along the electric held: 

F = Fn, = 1. (3.3) 


9 








In coordinates, where the unit vector n is the A^-th basic vector, the perpendicular components (3.2) acquire the form 
of the SO{N — 1) algebra standard basis. Among them, one can choose the — 1 commuting constants of motion. 

Apart from the angular momentum, the Coulomb-Stark system possesses a constant of motion, which is inherited 
from the TV-dimensional Runge-Lenz vector (2.16). It reads (see Ref. 24 for the TV = 3 case) 

p 

A = n ■ A — — (r^ — [n ■ x)^) (3.4) 

and complements the Liouville integrals to the full set. 


B. Integrable Hamiltonian of the Coulomb-Calogero-Stark system 

The inclusion of the inverse-square Calogero potential breaks down the SO{N) rotational symmetry of the initial 
Hamiltonian (2.1) to the discrete group containing the coordinate permutations (2.8b). Therefore, the different field 
directions n become nonequivalent. Among them, there is a preferable direction, which ensures the permutation 
invariance of the Stark term F ■ x: 

n=^(l,...,l). (3.5) 

This is virtually the direction of the center-of-mass coordinate, so that in the Jacobi system (2.32a), the Stark 
interaction takes the following form: 


F ■ X = Fxo ■ 


(3.6) 


So the generalized Calogero-Coulomb Hamiltonian 

Tfs™ = TTf" + Fxo 


in an external electric field is set to 
^2 


TT 


_7 

2 r 
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Fxo 


P I ^ 9{9-s^J) 7 p 


i<j 


^3? 


(3.7) 


We see that the particular choice of the field direction (3.5) keeps the entire Hamiltonian invariant under the coordinate 
permutations: 


(3.8) 

Therefore, it is well defined on the symmetric wavefunctions, where it is reduced to the Hamiltonian, as mentioned 
in the Introduction [see Eq. (1.2)]. We rewrite it here in a more compact form: 




, V 9(9-1) 

2 


7 

- -I- Fxq. 


(3.9) 


C. Constants of motion for "H®®" 

It turns out that the symmetry generators (3.2) and (3.4) to the Stark Hamiltonian (3.1) are extended straightfor¬ 
wardly to the generalized Calogero-Coulomb-Stark problem. 

First, consider the invariants of the generalized Calogero Hamiltonian Substituting the highlighted direction 

(3.5) into the orthogonal angular momentum tensor (3.2), which describe the remaining SO{N — 1) symmetry, we 
immediately get 

L^j = Lij -(- — 'y^{Ljk — Lik). (3.10) 

k 

Of course, this definition remains valid for the Dunkl angular momentum (2.11) too. Below we will consider just this 
case. It is easy to verify the following commutation relations among the coordinates and Dunkl angular momenta: 

— tXjSil xxiSji^ [T^j, :ro] — 'lipi (^*^^) 
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They ensure that the orthogonal components of the Dunkl angular momentum are constants of motion of the gener¬ 
alized Calogero-Coulomb-Stark problem: 




^7,FJ ~ 


(3.12) 


Alternatively, one can switch to the relative Jacobi coordinates and use the basis formed by the Lij (2.42b) for the 
description of the same algebra (3.10). Using the established identity S^i = Soi it is easy to ensure an analog of the 
commutation relations (3.11): 

- lytSji, [Lij,yo] = 0. (3.13) 


The above relation implies 


[4-,7^s“]=0. (3.14) 

The deformed analog of the modified Runge-Lenz component along the field direction is a constant of motion, too, 

A = xl) (3.15) 

= Xo + ^) - {rpr + + 3^0) ■ ( 3 - 16 ) 

where Ai = according to the definition (2.26), and we have used the identity (2.30). 

The proof is more complicated, and we present here some details of the derivation. Using the definitions (3.7), 
(3.15), and the already proven fact that the deformed Runge-Lenz vector is a constant of motion of the unperturbed 
Hamiltonian [see (2.18)], we obtain: 


=-^[Ai,a:o]-b j[7r^ -Xp]. (3.17) 

First, we calculate the second commutator from the right side of the above equation. Using a simple algebra and the 
relations (2.23), we get 

= — 2z(7r ■ X + X ■ tt) = —Airpr — 27V, (3.18a) 

[tt^, Xq] = [po, Xp] = -4zxopo - 2. (3.18b) 


Next, applying the identities (2.30) and (3.18b), the first commutator in Eq. (3.17) can be simplified as follows: 


--^[Ai,xp] = xp[7r^,xp] -b ixqPq + irpr -b 


7V- 1 


= i{rpr - XoPo) + 


N-l 


(3.19) 


Finally, substituting the commutation relations (3.18a), (3.18b), and (3.19) into the equation (3.17), we get the desired 
conservation condition: 


[A,77®“]=0. (3.20) 

It is easy to check that the modified Runge-Lenz longitudinal invariant commutes with the deformed angular momen¬ 
tum restricted to the relative coordinates: 


[LipA]=0 with l<i,j < N -1. (3.21) 

So we have proved that in the presence of a constant uniform electric field, the generalized Calogero-Coulomb model 
(3.7) still remains an integrable system. The integrals of motion are deformations of those for the conventional 
Coulomb model with the Stark interaction by the Dunkl operators. 
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D. Constants of motion of 

The integrals of the pure Calogero-Coulomb system (3.9) obtained by the restriction to the symmetric wavefunctions, 
must be symmetric too. Since the longitude component of the Runge-Lenz vector (3.15) obeys this condition, it remains 
as a correct integral for this system, 


= (3.22) 

We should take symmetric expressions of the kinematical constants of motion too, as in the absence of the electric 
field (2.25). For this purpose, it is more suitable to use the angular momentum in Jacobi coordinates (2.42a). Thus, 
we have 


['H^,F,t2k\ = Q, C2k= Y. (3.23) 

1<2<_7<A^ —1 

The first member of this set is the square of the relative Dunkl angular momentum; 

£2 = 2X+S'(S'-iV + 2). (3.24) 

Here I is the angular part of the Calogero Hamiltonian with reduced center of mass (the relative angular Calogero 
Hamiltonian). It also forms a single Casimir element of the relative Dunkl angular momentum algebra [9]: 

N — 2 X 

= 0, Ho = ^- i^—Py + (3.25) 

L L y 

Thus, we have proved the integrability of the Calogero-Coulomb-Stark system. 


E. Separation of variables in parabolic coordinates 

It is well known that the Coulomb-Stark system (3.1) admits separation of variables in parabolic coordinates, in 
which the Stark effect can be immediately calculated. Is this true for the Calogero-Coulomb-Stark system? Below we 
will show that the system admits complete separation of variables in parabolic coordinates for N = 2,3 and partial 
separation for N > 3. 

In the Jacobi coordinates (2.32a), the Calogero-Coulomb-Stark system acquires the following form, 

ny,F = 4- ^^— + Fxo+no. (3.26) 

^ v^o + y 

where again the last term is the Calogero Hamiltonian with the reduced center of mass (2.45). From the spherical 
coordinates {r,0,(pi,) (2.49), we move to the parabolic coordinates (C, ? 7 , using the map which leaves the relative 
angular variables ipi unaltered: 


^ = r-|-a;o, rj = r — Xq with Xo=rcos9. 
The inverse transformation reads 

a:o = ^—, = y = 

In the new coordinates the Hamiltonian (3.26) acquires the following form, 

il + + Br,) + — + -(p - rj), 

where 


_ 1 d d 


(3.27) 


(3.28) 


(3.29) 


(3.30) 
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We further proceed by extending straightforwardly the steps, applied for the usual Coulomb system in an external 
field in Ref. 25. Employing the following ansatz to the total wavefunction, 

V, ‘fi) = (3-31) 

we decouple the Schrodinger equation into three parts. Two of them depend on ^ and 77 , respectively, 

and replace the aforementioned equations (2.53a) and (2.53b) for the Calogero-Coulomb system in the spherical 
coordinates, 

0, (3.32a) 

0. (3.32b) 

Here the Coulomb charge is fractioned into two parts, 

Ai + A 2 = y. (3.32c) 


E 


^^+ 2 ^- 4 ^ - 


g(g + 7V-3) 


4^ 


+ Ai ) 4)i(^) = 


E E 2 9(9 + ^ —3) \ * / N 

+ - ^ -+ A2 ) Mv) = 


The third equation determines the eigenstate of the algular Hamiltonian of the relative Calogero model with the 
eigenvalue Iq given by (2.54). It remains unchanged and is given by Eq. (2.53c). 

Evidently, in Eqs. (3.32a) and (3.32b) the entire state (3.31) can be used instead of the partial wavefunctions. 
So, substituting $12 —t 4', one can eliminate E by subtracting from the first equation multiplied by 77, and from 
the second one multiplied by This manipulation together with Eq. (3.32c) yields the eigenstate equation for the 
modified Runge-Lenz invariant for the Calogero-Coulomb-Startk problem (3.15) in parabolic coordinates. 


H4' = (A2-Ai)4', (3.33) 

Therefore, the invariant A is responsible for the separation of the variables ^ and 77 in the parabolic coordinate system. 
It is used instead of the angular Hamiltonian I, which separates the radial and azimuthal variables in the spherical 
coordinates. The second invariant, given by the relative angular Hamiltonian I, is common in both cases and separates 
the relative angular degrees of freedom. 

In the absence of the electric held, F = 0, the fractional Coulomb charges Ai with 7 = 1,2 take discrete values 
depending on the parabolic quantum numbers rii = 0 , 1 , 2 ,...: 


_ rii + 5(9 + 1) 
"*'n + i(7V-3) 


(3.35) 


They also specify the partial wavefunctions 4)i(^) = $„^,j(^) and 4 * 2 ( 77 ) = ^n 2 qi'n)- Ttie principal quantum number, 
which characterizes the spectrum (2.56), is expressed now via parabolic quantum numbers: 


71 = 771+772+9+1- 


(3.36) 


As in the usual Coulomb problem [ 1 ], the electric held completely removes the degeneracy in the orbital momentum 
but preserves the degeneracy with respect to 9 . Using the similarity between the Calogero-Coulomb-Stark and TV¬ 
dimensional Coulomb-Stark problems, we can immediately write down the hrst-order energy correction caused by the 
weak electric held E in perturbation theory: 

^ + —y-) - ” 2 )F . (3.37) 


Hence, as in the spherical case, inclusion of the Calogero term preserves the perturbative spectrum of the system, but 
changes its degree of degeneracy. 


IV. TWO-CENTER CALOGERO-COULOMB SYSTEM 
A. Two-center Coulomb system 

As is known, the two-center Coulomb problem is integrable in any dimension [13, 26]. Let us locate the two central 
Coulomb charges 71^2 at the points a;i = a and X 2 = —a. The distances from them to the given space point are 
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denoted by ri and r 2 , respectively. Then the Hamiltonian acquires the following form: 


■lyo 7i 72 7i 72 . ^ 

2 |a:-a| |£c + a| 2 n rz' ^ 

The Hamiltonian (4.1) possesses the same spacial symmetry as the Stark problem described by the Hamiltonian (3.1). 
It remains invariant under the SO{N — 1) rotations, which preserve the line containing both charges. Let us denote 
a unit vector along this line by n: 

a = an, n? = 1. (4.2) 

According to the aforementioned symmetry, the orthogonal components of the angular momentum tensor L^-, given 
by (3.2), are preserved. They provide the system with A" — 1 commuting integrals of motion. 

The remaining A^-th constant of motion, like in the case of the Stark system, is provided by the suitably modified 
longitudinal component of the Runge-Lenz vector [13]: 

A = L'^ + (a ■ pY — 2(a ■ x) (— — —^ . (4.3) 

V ’’1 ^2 / 

In three dimensions, a close expression for this integral had been obtained before [27]. (See Ref. 28 for relations 
between both integrals.) 

It is easy to see that the above quantity commutes with . 


B. Two-center Coulomb system with Calogero term 

Consider now this system in the presence of the inverse-square Calogero potential. In order to construct the 
Hamiltonian of this system, we should replace the momenta operators by the Dunkl momenta as in previous sections, 
and then restrict the Hamiltonian to the symmetric wavefunction. In order to provide the system with permutation 
symmetry, like for the Stark Hamiltonian, we choose the line, connecting two Coulomb charges, to be directed along 
the center-of-mass coordinate (3.5). Hence, in the Jacobi coordinates (2.32a), the distances to the charges are given 
by 


n = \x - a\ = + (xq - a)^, r 2 = \x + a\ = ^Jy'^ + (xq -f a)^. 


(4.4) 


The generalized two-center Calogero-Coulomb Hamiltonian is 


njgen 

7i»72 


7l 72 
2 ri r 2 


(4.5) 


The particular choice of the direction of the vector a guarantees its permutation invariance. On the symmetric 
wavefunctions it gives rise to the two-center Calogero-Coulomb system defined in the Introduction (1.3), which we 
present in a compact form, 


= I g(g ~ 1) 

2 ^0? ri r2 ■ 


(4.6) 


Like the Calogero-Coulomb-Stark Hamiltonian, it possesses the symmetry given by the deformed angular momentum 
generators perpendicular to the predefined direction (3.10). 

It turns out that the modified Runge-Lenz integral (4.3) of the g = 0 Hamiltonian (4.1) is extended straightforwardly 
to the current case: 


A = £2 + o^pI 




(4.7) 


The proof is based on the extension of the elegant trick performed in Ref. 13, where it is shown that the Hamiltonian 
(4.1) preserves the modified Runge-Lenz component (4.3). 

First, we define the generalized Calogero-Coulomb Hamiltonians (2.7) with the shifted central charges: 


■ygen _ ^ _7l 

2 n ^ 


p/gen _ ^ _ 

2 r 2 ■ 


(4.8) 
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Then we combine in a proper way the symmetries of both Hamiltonians, which we have already discussed, in order 
to establish the conservation of the element (4.7) for the Hamiltonian (4.6). 

Each Hamiltonian (4.8) preserves the deformed angular generators, defined with respect to their Coulomb centers. 
Note that the Dunkl operators (2.8a) remain unchanged under the coordinate shift a; —>■ a; =F a. In turn, the deformed 
angular momenta are transformed under the shift of the original point according to the standard rule. For the first 
Hamiltonian, it reads as 




y/N 


TT,- — a:, — 


y/N 


T^i — Lij + (iTi TTj ) . 


(4.9) 


This provider the transformation rule for the deformed angular momentum square £ 2 , defined in (2.25): 


^ 2 ' = =£2 + - ’"j) 


i<j 


a 

m ■ 


'■,3 


= £2- 


- £0) ■ 


(4.10) 


The second term enters also in the expression of the longitudinal component of the deformed Runge-Lenz vector A\ 
(2.30). It behaves under the coordinate shift in the following way: 


Y - \/Na{Tr^-pi). (4.11) 

Hence, the entire component behaves under the coordinate shift as 

Xq — a 


- '/Nain'^ -pI) - //Vyi- 




ri 


(4.12) 


Take now a combination of these two invariants of the Hamiltonian H®™, in which the term with the anticommutator 
is eliminated: 


2 a 


Xq — a 


Cl^ + = ^2 - a^(7r^ - pI) - 2aji 

Substituting the expressions from (4.8), (4.10), and (4.12) into the conservation condition 


, £]^ + ^AT 


= 0 , 


(4.13) 


(4.14) 


and canceling out the terms containing the on the right part of the commutator, we come to the equation 

^2 


7r‘ 7 i ^ , 2 2 2071X0 

- -, £2 + a po- 

2 ri ri 


= 0 . 


(4.15) 


Since the square of the Dunkl momentum commutes with the deformed angular momentum (2.14), we have: 

[7r2,£2]=0. (4.16) 

This simplifies the commutator in the left side of the relation (4.15), and we can replace it by 


TT^ 2071 X 0 

+ 

7l ^ 1 2 2 

- , £2 + 0 Po 

[2 Xl J 


Ld j 


= 0 . 


(4.17) 


Remember now that the two Hamiltonians in Eqs. (4.8) are distinct from each other by the sign of the shift parameter 
and the values of the Coulomb coupling: 


O/gen _ O/gen I 
9,72 9.71 I 


>-a, 7i->72 ■ 


Therefore, an analogue of the equation (4.17) for 77®™ would be 



2072 X 0 


X 2 


+ 


72 ^ 2 2 

- , £2 + O Po 

X 2 


= 0 . 


(4.18) 


(4.19) 
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The sum of both equations produces the relation 


TT" / 71 72 

— , 20X0- 

2 V ?’! ^2 


7l , 72 ^ , 2 2 

-1-, £2 + a Pq 

n r2 


= 0 . 


(4.20) 


Again using the relation (4.16), we conclude that the above equation is equivalent to the conservation condition for 
the modified Runge-Lenz vector’s component along the symmetry axis of the system (4.7): 

,, A] = 0. (4.21) 


C. Separation of variables in elliptic coordinates 

Now let us show that in complete analogy with the previous case, the two-center Calogero-Coulomb system (4.6) 
admits complete separation of variables in the elliptic coordinates for A" = 2,3 and partial separation for N > 3. 

The map from the Jacobi variables {xo,y) to the elliptic coordinates (^, ?;) looks as follows (similar to the usual 
hydrogen atom case [25, 29]): 




ri +r2 
2a 


V = 


n - r2 
2a ’ 


(4.22) 


where is the distance from the i-th Coulomb charge (4.4). The relative angles pi remain unchanged. The new 
coordinates belong to the regions > 1 and — 1 < 7 < 1. The inverse transformation reads 


The two-center Calogero-Coulomb Hamiltonian (4.6) in elliptic coordinates reads 




1 


2a2(f2_^2) 


{Brj - B^) -I- 




7i 


72 


a{^ + v) 


where the operator B^ from the kinetic energy part acquires the following form: 

B, = _i -1 )^- 

Then, choosing the wavefunction, 

we can (partially) separate the variables in the Schrodinger equation = £’4' into the following parts: 

n(n 4- N - .71 

Be- 


(4.23) 


(4.24) 


(4.25) 


(4.26) 


£jj - 


~ ~^2 1 + 2a(7i-k 72 )^ + A^$i( 0 - 0 , 

(4.27a) 

~ ~^2 X +2a(7i ^ 2)11 +‘2-a?Eri^ A^ 4>2(?7) - 0 , 

(4.27b) 

(^ 2 i(v 3 ,, -q{q + N - 3)^ = 0- 

(4.27c) 


The third equation, which appears also in the spherical and parabolic cases (2.53c), describes the energy eigenstates 
of the relative angular Calogero Hamiltonian and its spectrum, depending on the composite quantum number q (2.54). 
In the absence of the Calogero term, it determines the spectrum and energy states of a free particle system on a 
{N — 2)-dimensional sphere. 

Obviously, the partial states <i>i 2 in the first two equations depend on the energy level E and the q. The parameter 
A in the first two equations separates the variables ^ and rj. It coincides with the eigenvalue of the slightly redefined 
Runge-Lenz invariant for the two-center Calogero-Coulomb system (4.7) with the Dunkl angular momentum square 
replaced by the doubled angular Calogero Hamiltonian, 

A = 21 + o^pI - 2axo ( ^ ^ V (4.28) 

V ?'i r2 J 
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The redefinition neglects the nonessential permutational invariant S term in Eq. (2.27). In the g = 0 limit, both 
integrals (4.7) and (4.28) become identical. In elliptic coordinates, the above formula looks like 


A = 


e-n 


2 


+ 


2{e + ri^-l) 
(^2 _!)(!_ ^ 2 ) 






(4.29) 


We follow the steps made above for the parabolic case. First, use the total wavefunction dt instead of the partial 
ones, ‘I’i_ 2 , in Eqs. (4.27a) and (4.27b). Next, cancel out the energy E by taking appropriate combinations of both 
equations. This yields the eigenstate equation for the modified Runge-Lenz invariant (4.29), 


= A4'. 


(4.30) 


It separates the two equations (4.27a) and (4.27b) from each other, as in the absence of the Calogero potential [25]. 


V. CONCLUSION 

In the present paper, we suggested new integrable generalizations of the rational Calogero model extended by i) the 
Coulomb potential with the Stark term and ii) the two-center Coulomb potential. We also found explicit expressions 
of complete sets of their constants of motion formulated in terms of the Dunkl operators. We demonstrated that these 
systems admit partial separation of variables in parabolic and elliptic coordinates, similar to conventional Coulomb- 
Stark and two-center Coulomb problems. This allows us to conclude that the suggested Calogero extensions preserve 
the spectra of initial systems but change the degree of degeneracy. As an illustration, we have written the first-order 
energy correction of the spectrum of then Calogero-Coulomb Hamiltonian caused by an electric field. 

The physical significance of the Coulomb-Stark and two-center Coulomb systems has been well known since (at 
least) the 19th century. The importance of the Calogero model does not raise any doubts either. The integrable 
synthesis, proposed in our paper, significantly extends its area of application. The study of some related systems is 
still in order. First of all, it is unclear whether similar integrable systems associated with an arbitrary root exist and, 
if so, how their hidden constants of motion look. But even in that case, we have no doubts that such generalizations 
will not admit the partial separation of variables. We intend to investigate these problems in a separate work. 
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